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ABSTRACT 

Red Supergiants (RSGs) are cool (~ 4000K), highly luminous stars (L ~ 10 5 L©), and are among the 
brightest near-infrared (NIR) sources in star-forming galaxies. This makes them powerful probes of 
the properties of their host galaxies, such as kinematics and chemical abundances. We have developed 
a technique whereby metallicities of RSGs may be extracted from a narrow spectral window around 
1 fim from only moderate resolution data. The method is therefore extremely efficient, allowing stars 
at large distances to be studied, and so has tremendous potential for extragalactic abundance work. 
Here, we present an abundance study of the Large and Small Magellanic Clouds (LMC and SMC 
respectively) using samples of 9-10 RSGs in each. We find average abundances for the two galaxies of 
[Z]lmc = —0.37 ±0.14 and [ZJsmc = —0.53 ±0.16 (with respect to a Solar metallicity of Z@= 0.012). 
These values are consistent with other studies of young stars in these galaxies, and though our result for 
the SMC may appear high it is consistent with recent studies of hot stars which find 0.5-0.8dex below 
Solar. Our best-fit temperatures are on the whole consistent with those from fits to the optical-infrared 
spectral energy distributions, which is remarkable considering the narrow spectral range being studied. 
Combined with our recent study of RSGs in the Galactic cluster Per OBI, these results indicate that 
this technique performs well over a range of metallicities, paving the way for forthcoming studies of 
more distant galaxies beyond the Local Group. 


1. INTRODUCTION 

The observed relationship between a galaxy’s stel¬ 
lar mass and its central metallicity, as well as the 
abundance trends as a function of galactocentric dis¬ 
tance in spiral galaxies, have provided vital clues as 
to how gal axies form and evolv e both in the l ocal uni¬ 
verse (e.g. iZaritskv et al.l 11994 iGarnettl I2002H and at 
larger redshift ( e.g. Tremonti et al.ll2004 lErb et al.ll2005 
iMaiolino et al.l I2008H . These observations have been 
used to test the theoretical predictions of various as¬ 
pects of galaxy formation and evolution under the frame¬ 
work of a dark energy and cold dark matter dominated 
universe, such as hierarchical clustering, infall, galactic 
winds and v ariations in the stellar Initial Mass Func¬ 
tion (IMFl dDe L ucia et all 12004 iKonnen et al.l 120071 : 


iColavitti et al.l 120081 : [Dave et al.l 120111 plus many oth¬ 
ers). 

However, deriving a galaxy’s abundances is a non¬ 
trivial task. Most commonly it has been done by measur¬ 
ing the strengths of Hil-region emission lines. A prob¬ 
lem is encountered when metallicities approach the Solar 
value, in that the temperature-sensitive ‘auroral’ lines 
become extremely faint, and instead one must rely on 
empirically calibrated ratios of strong lines. However, 
many such calibrations exist, and they are known to have 
large systematic differences in particular at high metal¬ 
licities. This has a profound effect on the reliability of the 


mass-metallicity relation and internal abundance gradi¬ 
ents of galaxies, both crucial diagnostics of galaxy evo¬ 
lution dKewlev fc Ellisonll2Q08l [Bresolin et al. [2009h . 

A promising alternative metallicity probe to H II- 
regions are evolved massive stars. They are highly lu¬ 
minous (~ 10 5 L©), allowing them to be studied at Mpc 
distances. Correct interpretation of their spectra is made 
possible by sophisticated model atmospheres, and so such 
work is free of the systematic effects that hamper stud¬ 
ies of Hll-regions such as unknown temperature struc¬ 
tures and large-scale inhomogeneities. A small sample of 
nearby gal axies have been studied using blue supergiants 
(BSGs, e.g. lKudritzki et al]l2012L [201311201^ ). and it was 
shown in the case of NGC 300 that the BSG results were 
in excellent agreement w ith those of aur oral line Hll- 
region measurements dBresolin et al.l[2009h . 

Another promising probe of extra-galactic abundances 
are Red Supergiants (RSGs). These stars are luminous 
{L .> lO 5 ^©), and with fluxes which peak at ~l/im are 
among the brightest NIR sources in a galaxy. Their 
brightnesses and colours make them extremely easy to 
identify, with several thousand such stars expected in a 
Milky-Way type galaxy. Their young ages (<; 20Myr) 
mean that their surface abundances are representative 
of those in the interstellar medium, aside from those of 
carbon and nitr ogen w hich are s lightly modulated by 
CNO burning (jDavies et all l2009h . Historically, abun- 
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dance studies of RSGs have required spectra with high 
resolving powers (A/AA = R ^ 20,000) in order to dis¬ 
tinguish the diagnostic metallic lines from the millions 
of overlapping mo lecular lines (e.g. iCunha et, al.l 120071 : 
iDavies et a l. 2009). This has limited their usefulness in 
extragalactic science, since prohibitively large integration 
times are requi red at distan c es bey ond IMpc. 

However, in [Davies et al. (20.10. hereafter Paper I), 
we demonstrated that by isolating a spectral window in 
the J-band which is relatively free of molecular lines, 
abundances may be extracted at resolutions as low as 
R ~3000. Not only does this dramatically reduce the 
required exposure times, but is also suitable for multi¬ 
object spectrographs. This then would allow entire 
galaxies at several Mpc to be abundance-mapped in only 
a few hours with an 8 m-class telescope. In a further 
paper, we showed that this technique was ideally suited 
to extremely large telescopes (ELTs), with tremendous 
gains due to the large aperture and the adaptive optics 
system optimised for the NIR. We estimated that with 
the E-ELT we could expect to obtain abundances from 
individual RSGs out to dist ances o f tens of Mpc in less 
that one night’s integra tion (|Evans et all l20lTh . 

In iGazak et al.l (|2014ll , we performed a comprehensive 
test of this technique within a Solar metallicity environ¬ 
ment, analysing a sample of stars in the Milky Way 
open cluster Per OBI. We found an average metal¬ 
licity Z = —0.04 ± 0.08, consistent with the results 
from blue supergiants within the_same clu ster a nd in 
the Solar neighbourhoo od dFirnstein fe P rzvbilla 2012; 
iNieva fe Przvbillal f2012l. We also demonstrated defini¬ 
tively that the technique works down to resolving powers 
less than R = 3000. In this paper, we undertake a similar 
analysis in two lower metallicity environments, the Small 
and Large Magellanic Clouds (SMC and LMC respec¬ 
tively), to test this technique for observations spanning 
a broad range of metallicities. 

We begin in Sect. [2] with a description of the data and 
its reduction. We describe the model grid we use, along 
with the analysis technique, in Sect. [3] The results are 
presented in Sect. [2 and we compare our results to other 
studies of massive stars in Sect. [5] We conclude in Sect. 

I 

2. OBSERVATIONS AND DATA REDUCTION 
2 . 1 . Observations 

We obtained observations of several s tars in the L MC 
and SMC using VLT+XSHOOTER ()D’Odorico et all 
120061 1 under ESO programme number 088.B-0014(A) (PI 
B. Davies). These observations provide continual spec¬ 
tral coverage from 0.3-2.4/mr. The description of the ob- 
servin g strategy and the sample selection was described 
in IDavies et all (120131 , hereafter D13), which we briefly 
summarise here. The stars were observed in nodding 
ABBA mode with at least four exposures per star, and a 
randomized jitter at each position on the slit. For each of 
the instrument arms - UVB, VIS, and NIR - we used the 
5.0" slits to minimize slit losses and obtain accurate spec- 
trophotometrj0. The use of the broad slit meant that 
spectral resolution was determined by the seeing, which 

1 The spectrophotometric aspect of these observations was not 
required in the present study, but was essential for the study pre¬ 
sented in D13. 


at 1" corresponds to a resolution R = AA/A ~ 5000. 
The precise value of R for each star was determined at 
the analysis stage (see Sect. 0 , and was found to be 
within R = 4000 — 8000. Integration times were chosen 
to achieve a signal-to-noise ratio (SNR) of at least 100 
in the J-band (see Table 1 of D13). Flux standard stars 
were also observed each night and to correct for the at¬ 
mospheric absorption in the NIR, telluric standard stars 
of spectral type late-B were observed within one hour of 
each science target. In general observing conditions were 
good on each night, and the seeing was around 1 " or 
better. The standard suite of XSHOOTER calibration 
frames used in the data reduction process were taken 
at the beginning and end of the night (for details see 
IModigliani et alJl2010D . 

2 . 2 . Data reduction 

The initial steps of the reduction process were 
done using the XSHO OTER data reduction pipeline 
(IModigliani et al.l [20104 . These steps included subtrac¬ 
tion of bias and dark frames, flat-fielding, order extrac¬ 
tion and rectification, flux and wavelength calibration. 
The accuracy of the wavelength solution was determined 
by measuring the residuals of the fitted arc lines with 
their vacuum wavelengths. 

The spectra of the science targets and the telluric 
standards were then extracted from final rectified two- 
dimensional orders. The strengths of the telluric lines 
were scaled in order to give the best cancellation across 
the diagnostic bandwidth (1.15-1.22/rm). 

3. ANALYSIS 
3.1. Model grid 

For our project we have computed a new grid of model 
atmospheres. The se atmospheres were ge nerated using 
the MARCS code (iGustafsson et all 12008 1. which oper¬ 
ates under the assumptions of LTE and spherical sym¬ 
metry. The grid of models is four-dimensional, computed 
with a range of metallicities (log Z), gravities (log < 7 ), ef¬ 
fective temperatures (T e fj) and microturbulent velocities 
(£). The chemical composition was scaled from Solar 
between [Fe/H]=-1.5 and +1.0 in steps of 0.25dex; T e g 
between 3400 and 4000K in steps of 100K, with fur¬ 
ther models at 4200K and 4400K; logg between —1.0 
and +1.0 in steps of 0.5 (in cgs units); and values of £ 
from 1 to 6 km s . All models were computed with an 
adopted stellar mass of M*=15M 0 . Though RSGs may 
have masses between ~8-25M 0 , the only effect on the 
model of changing M* is to alter the pressure scale height, 
which remains largely unchanged throughout the RSG 
mass range (see discussion in Paper I). T he S olar-scaled 
abundance ratios were taken from lGrevesse et all ( 2007 1. 
The synthetic spectra were computed using the updated 
versio n of the SIU code, as described in lBergemann et all 
(I2012H . SIU includes new opacities which are consis¬ 
ten t with those in the DETAIL code dButler fe Giddingsl 
H985i l which we use for statistical equilibrium calcula¬ 
tions of Ti, Mg, Fe, and Si. 

To increase the precision of our abundance measure¬ 
ments we have undertaken the task of computing correc¬ 
tions to our diagnostic lines which take into account de¬ 
partures from local thermodynamic equilibrium (LTE). 
Corrections to the Fel, Til and Sil lines are presented 
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in iBergemann et all (12012L I2013D . Corrections to the 
Mg I lines are still in preparation, and are not consid¬ 
ered in this work. The list of diagnostic lines used in 
this study are listed in Table HI For this current study 
we have avoided using the Fei lines at 1.163826/im and 
1.15936/zm as these lines were subject to strong contami¬ 
nation from telluric absorption, making the line strengths 
uncertain. 


TABLE 1 

Diagnostic spectral lines used in this analysis 


Species 

A(p,m) 

E\ow{eV) 

■Eup(eV) 

log (s/) 

Fei 

1.188285 

2.20 

3.24 

-1.668 

Fe I 

1.197305 

2.18 

3.22 

-1.483 

Si I 

1.198419 

4.93 

5.96 

0.239 

Si i 

1.199157 

4.92 

5.95 

-0.109 

Si i 

1.203151 

4.95 

5.98 

0.477 

Si i 

1.210353 

4.93 

5.95 

-0.351 

Til 

1.189289 

1.43 

2.47 

-1.730 

Til 

1.194954 

1.44 

2.47 

-1.570 


3.2. Continuum placement 

Our method is based on matching the observed line 
strengths relative to what we deem to be the conti n- 
uum level. As described in Paper I and in lEvans et al.1 
(12011 . accurate continuum placement for RSGs is non¬ 
trivial. Their optical and near-IR spectra contain many 
thousands of blended molecular absorption lines, mean¬ 
ing that the true continuum may be well above the lo¬ 
cal maximum observed flux level. The spectral win¬ 
dow that we adopt for this work was chosen specifically 
because of the low contamination by these lines, mak¬ 
ing continuum placement easier. Nevertheless there is 
still some weak molecular absorption present, and the 
‘pseudo-continuum’ we observe may be a few percent be¬ 
low than the true continuum. However, as long as the 
methodology we employ to determine the level of the 
pseudo-continuum is the same for both models and data, 
we will always find th e same best-fitt ing model (for fur¬ 
ther discussion, see lGazak et all 12014 1. 

The first step of this process is to flatten out any large- 
scale slope that the data may have within the J-band 
window. To measure and remove this slope we first di¬ 
vide through by a template model spectrum which has 
had the strong absorption lines masked out. We use a low 
metallicity template, since these spectra have the clean¬ 
est continua free of molecular absorptiorQ. The resulting 
ratio spectrum was heavily smoothed with a low-order 
median filter, before fitting a 2nd-order polynomial. This 
polynomial is our measure of the low-order ‘tilt’ of the 
observed spectrum, which is then removed by dividing 
through by the polynomial. 

To find the continuum level of this flattened spectrum, 
we ranked the pixels in order of flux values and deter¬ 
mined the continuum to be the median of those with 
the top 20%. This exact same methodology was applied 
to both data and models, to ensure that the placement 

2 We investigated the effect on our results of using different tem¬ 
plates, and found there was no effect p rovide d the template had 
metallicity lower than [Z]=0.0, see Seel . IS.4.41 


of the pseudo-continuum was consistent. We investi¬ 
gated the robustness of our results to the exact method 
of continuum placement, experimenting with (for exam¬ 
ple) varying the cutoff value for selecting pixels to fit, 
varying the parameters of the template spectrum in the 
tilt-correction step (see above), as well as experiment¬ 
ing with alternative meth ods of continuum fitting such 
as those presented in lEvans et al.l (1201 il l. We found that 
our results are not sensitive to this, provided the same 
method is applied to model and data, and the pixel cut¬ 
off is greater than ~60% (i.e. we use the top 40% of 
ranked pixels) when determining the median continuum 
flux level. 

We also tested the sensitivity of our results to the 
signal-to-noise (S/N) of our data, by adding various lev¬ 
els of Gaussian noise. Since our method for continuum 
placement selects the highest value pixels, in the event of 
low S/N we would be selecting not the continuum but the 
high flux tail of the Gaussian noise. This would result in 
the continuum being placed too high, leading to the line 
strengths relative to the continuum being over-predicted, 
and ultimately an overestimate of the metallicity. We did 
indeed find that at low S/N the average abundance lev¬ 
els increased, but provided the S/N per pixel is greater 
than ~60 then the measured abundances are stable at 
the level of ±0.02dex. In practice, we have shown that 
it is possible to measure accurate abundances below this 
S/N limit by degrading the spectra to lower resolution 
at constant sampling of the resolution element, provided 
the degraded resolution is above R = 3 000 and the S/N 
per resolution element is ;> 100 dGazak et al.ll2014lf . 

3.3. Spectral fitting methodology 

The basis for our analysis is that described in Paper I. 
Briefly, we compare the strengths of absorption lines in 
an input spectrum with those of a template spectrum at 
each point in the model grid. The difference between the 
input and model spectra is used to compute a y 2 value, 
and the best fitting model is deemed to be that with the 
lowest y 2 . 

Before computing the y 2 value at each point in the 
grid, there are two subtle effects that must be taken into 
account which may produce spuriously high y 2 values. 
The first is that any relative velocity shift between model 
and data must be removed. This is done by an iterative 
cross-correlation procedure, with the model shifted until 
the measured offset is below 0.1 pixels. 

The second effect is that of a mismatch between the 
line-widths of the models and the data. Such a mis¬ 
match may be caused by variations in instrumental spec¬ 
tral resolution, due to e.g. instrument flexures or see¬ 
ing variations, or by astrophysical effects such as macro¬ 
turbulence and stellar rotation. Since the data pre¬ 
sented here were taken with a wide 5" slit, seeing vari¬ 
ations can cause the spectral resolution to vary from 
4000 < R < 8000. In an equivalent width (EW) study 
we would not need to account for such factors as the flux 
within each line would be conserved. However, in the 
case of RSGs, broadening the spectral lines means that 
more of the unresolved molecular absorption that con¬ 
tributes to the pseudo-continuum is swallowed up by the 
diagnostic lines, meaning that the measured equivalent 
width actually increases as the resolution is degraded. 

To account for this, we again employ an iterative pro- 
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cedure. First, we take a grid of model spectra that have 
been degraded to a resolution of R = 10,000, which we 
expect to be substantially higher than that of our data. 
We then find the model in this grid which has the best 
overall match to the diagnostic line strengths, i.e. the 
contrast between the the line-centre and the continuum. 
This quantity is obviously dependent on the spectral res¬ 
olution, as for a given equivalent width this contrast will 
be greater for spectra of higher resolving power. We 
then measure the widths of the diagnostic lines in both 
the best-fit model and the observed data by fitting Gaus¬ 
sian profile^. In this first iteration, we expect the line 
widths of the model to be narrower than the data, as we 
have deliberately overestimated the spectral resolution. 
We compare the difference in the fitted profile widths of 
the model and observed spectra, use these to estimate 
the true resolution of the data. This process is then re¬ 
peated with a model grid that has been degraded to the 
updated estimate of the resolution, and we continue to 
iterate until the input and output resolution is stable. 
Applying this methodology we find that the estimated 
resolution converges to 5% within four iterations. This 
corresponds to a systematic error in the inferred value of 
[Z] of ±0.04dex, which we account for in our uncertain¬ 
ties (see Sect. 13 43 . 

To check the sensitivity of our results to spectral res¬ 
olution, we took the observed data and degraded it to 
lower resolutions, and re-performed the abundance anal¬ 
ysis. We found that the best-fit parameters were sta¬ 
ble down to measured resolutions of R > 2000, at which 
point the diagnostic lines begin to blend together. This is 
agai n consistent w ith our results for the RSGs in Per OBI 
(iGazak et al.ll20lih 

3.4. Sensitivities of the diagnostic lines to the free 
parameters 

Since there are four model parameters which can po¬ 
tentially influence the output spectrum, it is important 
to understand how each of these variables affects the rel¬ 
ative strengths of the diagnostic lines we employ. Be¬ 
low we discuss the sensitivities of the diagnostic lines to 
each of the free parameters in our model grid. We sup¬ 
ply demonstrative figures (Figs. 1 - 3) to illustrate how 
certain lines respond to changes in one input parame¬ 
ter while the others are fixed, however we note that our 
fitting methodology is to consider all lines and all param¬ 
eters simultaneously (see Sect. 13.4.41) . 

3.4.1. Microturbulence, £ 

The effect of this parameter is to increase the EW of 
stronger saturated lines, whilst the EW of weaker lines 
will be less affected. We therefore expect that the rela¬ 
tive strengths of lines of the same element and ionization 
stage will be sensitive to this parameter. This can be 
simply demonstrated by studying the variation in line 
strengths as £ is increased while keeping the other pa¬ 
rameters fixed, where we see the two Fe I lines responding 
differently (Fig. [Q left panel). We note that the crossing 
of the two curves in the left part of Fig. 1 is a subtle 
non-LTE effect. Both lines belong to the same multiplet 

3 The dominant source of line broadening in our data is instru¬ 
mental, so Gaussian profiles are an appropriate choice for modelling 
the lines in our spectra. 


and, thus, in LTE the curves displaying their strengths 
would not cross. However, in our non-LTE calculations 
the multiplet transitions are treated as individual lines 
and have distinct departure coefficients which respond 
differently when the atmospheric parameters are varied. 

3.4.2. Effective temperatures T e g 

Spectra of RSGs in the J-bancl window contain lines 
that have different excitation potentials (Table [I] sec 
iBergemann et all 120121 120131 ). Therefore, we may ex¬ 
pect that the ratio of (for example) the Si I to Til line 
strengths would be sensitive to T e g. This can be demon¬ 
strated by looking at how the strengths of these lines 
change as a function of T e g with the other parameters 
fixed (Fig. [l] right panel). 

3.4.3. Gravity and metallicity 

Separating the effects of these two parameters is less 
straight-forward than for T e g and £. Overall, increas¬ 
ing (decreasing) metallicity (gravity) at fixed T e g and £ 
has the effect of increasing the strengths of all the lines, 
though some lines are more sensitive than others. This 
means that there is a degree of degeneracy between these 
two parameters. 

To investigate this further, we performed principle 
component analysis (PCA) on the matrix of the strengths 
of the eight diagnostic lines. We did this in two stages - 
first, we extracted two dimensional grids of line-strengths 
for each line as a function of T e ff and £, at fixed logg 
and [Z]. We computed the eigenvectors, and searched for 
pairs of eigenvectors which were orthogonal, preferring 
the earliest possible components since these account for 
the greatest variance. We then repeated this process al¬ 
lowing logg and [Z] to vary at fixed T e g- and f. 

We show an example of this in Fig. [2] The left panel 
shows the first pair of orthogonal eigenvectors when al¬ 
lowing T e s and f to vary. Here, the first two eigenvectors 
eO and el display a large degree of orthogonality, and 
which account for 99.8% of the variance. The interpre¬ 
tation of this is that these two parameters are readily 
constrained from the very different effect that each has 
on the strengths of the diagnostic lines. Overplotted on 
this panel is the datum for the first star in our sample, 
LMC 064048, for which we have computed the projec¬ 
tion onto the eO —el plane using the star’s line strengths. 
The very small error bars illustrate that T e s and £ can 
be constrained to high precision. 

In the right panel of Fig. [2] we show a similar anal¬ 
ysis but for varying logg and [Z]. Here it can be seen 
that to find orthogonal eigenvectors we had to go to sev¬ 
enth eigenvector e6, which accounts for only a very small 
fraction of the variance (<Cl%). This is illustrated by 
the error bars on the observation of LMC 064048, which 
are much larger relative to the parameter space covered 
by the models. The interpretation of this is that the de¬ 
generacy between log g and [Z] is harder to break due to 
the subtle differences that exist between their effects on 
the strengths of the spectral lines, and detecting these 
differences requires high signal-to-noise data. 

Since obtaining such high S/N data may be challeng¬ 
ing, especially when considering that there may be sys¬ 
tematic errors in both the data reduction and the spectral 
synthesis that we have not accounted for, we explored an 
alternative method to constrain logg and hence better 
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Fig. 1.— An illustration of the effects of varying the model parameters on the line strengths. In the left panel, we plot the line strength 
(= continuum - central line depth) for the strongest Fel lines as a function of £, demonstrating how the two lines behave differently. Right 
panel: the strengths of the Til and Sil lines as a function of T e ff. 



Fig. 2.— Principle component analysis of the line strengths (= central line depth) in our model grid, illustrating the orthogonality between 
pairs of model parameters. In the left panel, we plot the first two principle components (eO and el) when T e ff and £ are allowed to vary 
at fixed log g and [Z]. For demonstration purposes we have taken the observed line-strengths of the first star in our sample, LMC 064048, 
and used the eigenvectors to calculate eO and el (overplotted as a blue square). In the right panel we do the same, but for the other pair 
of model parameters. This time, the first pair of orthogonal components were eO and e6. 


understand the degeneracy between this parameter and 
[Z]. First, the model grid was subsampled in flux onto 
a finer grid, with 0.2 km s _1 steps in £, O.ldex steps in 
[Z], and 0.25dex steps in logg. The strengths of the di¬ 
agnostic lines were measured for each model in the grid. 
For each subsampled value of [Z], we extracted a three- 
dimensional sub-grid (T e ffj £, and logg), and performed 
a ^-minimization search to find the best fitting values of 
the three free parameters at that value of [Z]. Since the 
best-fitting T e ff and £ are well constrained for a given set 
of line-strengths (see above), this result then tells us the 


best-fitting value of log g for that particular value of [Z]. 
We then repeat this for all values of [Z]. 

An example of the results of this analysis can be found 
in Fig. [3l which shows the best-fit values of T e ff, logg 
and £ for a fixed value of [Z], again for the first star in 
our sample. The plot further helps to illustrate the de¬ 
generacies in our fitting method. The best-fitting values 
of the parameters T e g and £ are stable regardless of what 
value [Z] is fixed at, and hence that there is very little 
degeneracy between them. 

The second-top panel demonstrates the degeneracy be- 
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Fig. 3. — The best-fitting values of T e ff (top panel), log g (sec¬ 
ond panel), and £ (third panel) for fixed values of [Z] for the star 
LMC 064048. The bottom panel shows the x 2 at each [Z]. The dot¬ 
ted lines in the second panel illustrate the maximum and minimum 
possible values of log g , determined from the star’s photometry and 
temperature. The circled points are those selected for determining 
the x 2 -weighted average, since it is these points which are within 
the possible gravity range (or closest to it). The dotted lines in the 
top and third panels denote the average T e ff and £ respectively. 

tween [Z] and logg. We see that for low values of [Z], 
lower gravities are preferred. This is because lowering 
[Z] reduces the line strengths, and so to match the line 
strengths of the data we must select low values of log g 
to compensate. Similarly, at high metallicities, we must 
select the highest gravity models to best match the line 
strengths. At intermediate metallicities we see a transi¬ 
tion from low to high gravities. 

However, we can put prior constraints on logg. For 
a given T e g and luminosity (the latter obtained from 


D13) we can estimate the radius R of the star. Then, 
by selecting an appropriate mass range M for that 
value of L using evolu tionary models (8-30M 0 , e.g. 
iMevnet fo Maederl 12000( 1. we can estimate the surface 
gravity from g = GM/R 2 . Allowing for a contribu¬ 
tion from convective pressure which would redu ce the 
effective gravity by up to 0.3dex dChiavassa et all[20111 : 
lOazak et al.l 12014( 1, we can identify the extreme possi¬ 
ble values of logg. When taking the y 2 -weighted mean 
of each of the parameters within this range, or simply 
just their values at the y 2 minimum, we find very similar 
results to those in the by the PCA analysis (see above). 

3.4.4. Best-fit parameters and their uncertainties 

Our methodology for finding the best-fitting model is 
as follows. Using the subsampled grid (see above), we 
compare the line-strengths of the observed spectrum O 
with those in the model grid A4 at the best-fitting res¬ 
olution. At each point in the grid we compute the un¬ 
reduced y 2 value from the eight diagnostic lines, 

n 

X 2 = ( 1 ) 

i 

where a is the error on the line-strength, determined by 
the S/N, and which does not account for the error in the 
spectral resolution which is treated separately (see be¬ 
low). We ignore all points in the grid where the value 
of log g is unphysical. Since we know the luminosities of 
the stars in our sample, this also allows us to rule out 
grid points with low-T e fj / high-logg (and vice versa), 
since stars at these grid points would have anomalously 
high (or low) masses. The best-fitting values of each 
parameter are determined from the weighted mean of 
all points in the grid, where the weight w is given by 
Wi = exp(—y 2 ). Eliminating unphysical regions of pa¬ 
rameter space before computing the averages is useful 
for eliminating localised y 2 -minima, which may skew the 
best-fits for the parameters. 

When determining the uncertainties on each of the four 
free parameters, we consider three sources of random er¬ 
rors. We discuss each of these below individually. The 
total errors are taken to be the quadrature sum of these 
three sources of uncertainty. 

‘Degeneracy’ errors: — This is the term we use to describe 
the uncertainty caused by several points within the model 
grid having y 2 values which are close to that of the best¬ 
fitting model. To determine this error, we define a y 2 
tolerance, Ay 2 = Xmin + 3, where y 2 m i n is the minimum 
unreduced y 2 . The tolerance value of 3 is chosen for the 
following reasons. If all line strengths in the observed 
spectrum are on average matched to lex by those in the 
model spectrum, the unreduced y 2 value should be equal 
to the number of diagnostic lines n, which in this case 
n = 8. Should one of the lines be fitted to only 2cx, 
the rest being fitted to lex, the unreduced y 2 becomes 
(n — 1) x l 2 + 1 x 4 2 = n + 3 0. We consider all models 
which are able to fit the data to better than Ay 2 to define 
the limits of the degeneracy errors. 

4 This tolerance threshold for \ 2 compares well to the ‘classical’ 
value of A\ 2 =2.3 for a la deviation from the peak of a purely 
Gaussian probability distribution. 
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Fig. 4.— Planes through the model grid showing the degeneracy between the parameters for the star LMC 064048, see text for how 
planes were created. The black contours indicate contours of equal x 2 values, and are drawn at levels of Xmin ± (1) 2, 3, 5,10). The red 
contour is drawn at x 2= X 2 ■ +3, which indicates our la uncertainty. The dotted lines on the right-hand plot indicate the minimum and 
maximum possible values of'fog g given the star’s temperature and luminosity. 


To illustrate these errors, we calculate the y 2 values at 
every point in the grid and take two-dimensional projec¬ 
tions. For one pair of parameters, say T e s and [Z], we 
determine the minimum y 2 for every value of T e g and 
[Z] whilst allowing the other two parameters (in this case 
log<? and y) to be free, other than the limits we have al¬ 
ready placed on log g (see Sect. 13.4.TT1) . We thus construct 
a 2-D plane of minimum y 2 values as a function of T e g- 
and [Z]. On each of these planes, we then draw contours 
of equal y 2 . An example of this is shown in Fig. [4] This 
figure shows again how tightly constrained T e ff and £ are, 
whilst also demonstrating again the degeneracy between 
[Z] and log g. Typical errors on [Z] from this source of 
uncertainty are between 0.1 — 0.2dex. 

‘Resolution’ errors :— As discussed in Sect. ED9 part of 
the fitting involves determining the spectral resolution, 
which we do by matching the line widths at the observed 
line depths. We do this by applying an iterative pro¬ 
cedure until convergence at the 5% level is achieved. 
Were we basing our fitting methodology on line equiv¬ 
alent widths then this uncertainty would be combined 
with the degeneracy errors. However, since we are fit¬ 
ting both the line depths and the line widths, we consider 
these sources of error separately. 

To determine the magnitude of this uncertainty on our 
results, we re-ran the analysis above but compared the 
observed line strengths to those of models which had 
been degraded to the observed resolution ±5%. The ef¬ 
fect on T e ff and £ was minimal, since these parameters 
are sensitive to line depth ratios, whereas altering the 
resolution uniformly alters only the line depths. The ef¬ 
fect on log g was also small compared to the degeneracy 
errors, typically <^0.1dex. The effect on [Z] was again mi¬ 
nor, ±0.04 on average, but large enough that they cannot 
be neglected. 

Continuum placement errors: — These errors refer to the 
uncertainties described in Sect. 13.21 which concern the 
placement of the pseudo-continuum. As has already 
been discussed, this is not the true continuum since un¬ 
resolved molecular absorption reduces the overall maxi¬ 
mum flux level. However, as long as the same principles 


and methodology are applied to both models and obser¬ 
vations then the line to psuedo-continuum ratio is still 
an accu rate measu re of metallicity (see also discussion in 
iGazak et al.ll20lJ) . 

To investigate the effect of this uncertainty we re-ran 
the analysis several times, each time fitting the contin¬ 
uum in a slightly different way. We varied aspects such 
as the metallicity of the template spectrum (used to cor¬ 
rect for any spectral tilt), the fraction of pixels used to 
determine the highest flux level, as well as alternative 
continuum placement methods we have e xperimented 
with previously, such as those described in lEvans et al.l 
(120111) . Each time the difference in the fitted parame¬ 
ters was minimal, with the abundances stable to within 
±0.02dex. We conclude that this source of error is small 
compared to the other two described above, but is in¬ 
cluded nonetheless. 

4. RESULTS OF THE J-BAND ANALYSIS 

The data, along with the best fitting model spectra, are 
shown in Figs. [5] and [G] While the match to the diagnos¬ 
tic lines is excellent, the fits to the unresolved features of 
the pseudo-continuum is also very good. Since the latter 
were not used to diagnose the best fitting model param¬ 
eters, this gives us further validation that our models are 
performing well. The only features not well matched are 
the two Mg I lines, and an unknown feature at 1.2045/rm. 
The explanation for the former is likely non-LTE effects, 
the corrections for which will be presented in a forthcom¬ 
ing paper (Bergemann et ah, in prep), while we currently 
do not have an explanation for the latter. 

The best-fitting model parameters are listed in Table 
[21 As described in Sect. 13.41 the parameters T e ff and £ 
are sensitive to line ratios, and with high S/N data in 
which the lines are well resolved, these parameters are 
strongly constrained. As discussed in the previous sec¬ 
tions and illustrated in Figs. [3] and [4) T e g and £ are stable 
to well within ~10% regardless of what value we fix the 
metallicity to. The major source of error in the metallic¬ 
ity [Z\ is the degeneracy with log g, which accounts for 
an uncertainty of approximately ±0.1dex. 

For the LMC stars, we find an average metallicity of 
[Z\ = —0.37 ± 0.14, consistent with other measurements 
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TABLE 2 

Best fit model parameters for each star in this study. 


Star 

T e s (K) 

l°g(g) 

£ (km 

-/b) 

[Z] 


R{± 5%) 

LMC-064048 

3860 

± 

70 

-0.2 ± 

0.4 

3.1 

± 

0.2 

-0.42 

± 

0.17 

5600 

LMC-067982 

3910 

± 

60 

-0.3 ± 

0.3 

3.7 

± 

0.3 

-0.43 

± 

0.14 

5400 

LMC-116895 

3950 

± 

60 

-0.3 ± 

0.5 

3.2 

± 

0.2 

-0.30 

± 

0.18 

6400 

LMC-131735 

4110 

± 

50 

-0.0 ± 

0.2 

3.3 

± 

0.2 

-0.50 

± 

0.09 

6700 

LMC-136042 

3850 

± 

50 

0.1 ± 

0.1 

2.7 

± 

0.2 

-0.07 

± 

0.07 

6900 

LMC-137818 

3990 

± 

50 

0.0 ± 

0.2 

2.3 

± 

0.2 

-0.50 

± 

0.13 

8200 

LMC-142202 

4100 

± 

50 

-0.5 ± 

0.3 

4.2 

± 

0.2 

-0.28 

± 

0.08 

5100 

LMC-143877 

4060 

± 

80 

-0.3 ± 

0.5 

3.7 

± 

0.3 

-0.29 

± 

0.17 

4000 

LMC-158317 

4040 

± 

80 

-0.0 ± 

0.5 

3.6 

± 

0.4 

-0.37 

± 

0.19 

4600 

Average, LMC 









-0.37 

± 

0.14 


SMC-011709 

4020 

± 

50 

-0.3 ± 

0.1 

3.0 

± 

0.2 

-0.53 

± 

0.10 

4900 

SMC-013740 

3970 

± 

70 

0.3 ± 

0.5 

2.6 

± 

0.2 

-0.59 

± 

0.22 

5800 

SMC-020133 

3970 

± 

80 

-0.2 ± 

0.4 

2.3 

± 

0.3 

-0.21 

± 

0.20 

5000 

SMC-021362 

3970 

± 

70 

0.0 ± 

0.5 

2.7 

± 

0.2 

-0.57 

± 

0.20 

5400 

SMC-030616 

4040 

± 

60 

-0.2 ± 

0.5 

2.8 

± 

0.2 

-0.47 

± 

0.17 

7000 

SMC-034158 

4180 

± 

70 

0.1 ± 

0.5 

3.2 

± 

0.2 

-0.48 

± 

0.15 

7300 

SMC-035445 

4040 

± 

80 

-0.0 ± 

0.5 

2.4 

± 

0.2 

-0.53 

± 

0.21 

5200 

SMC-049478 

4110 

± 

60 

-0.3 ± 

0.4 

3.2 

± 

0.2 

-0.24 

± 

0.12 

4900 

SMC-050840 

3920 

dfe 

50 

-0.3 ± 

0.3 

3.0 

± 

0.2 

-0.72 

± 

0.12 

5400 

SMC-057386 

4120 

± 

50 

-0.0 ± 

0.5 

2.9 

± 

0.2 

-0.57 

± 

0.15 

5900 

Average, SMC 









-0.53 

± 

0.16 




Fig. 7.— Metallicity distributions of the RSGs in the two Mag¬ 
ellanic Clouds. 

of young stars in this galaxy (see following section). The 
one object which is slightly discrepant with the others is 
LMC 136042 which has a metallicity 0.3 dex higher than 
the average. As discussed in D13, this star’s spectrum 
is blended with a nearby hot star, evident from the blue 
region of the spectrum. This companion should not con¬ 
tribute much to the J-band flux, indeed if it did then we 
would expect the diagnostic lines to be diluted and for 
the fitted metallicity to be abnormally low. Therefore, 
we do not see a good reason to discard this object from 
our sample. 

In terms of the SMC stars, we find an average metal¬ 


licity of [Z] = —0.53 ± 0.16, slightly on the higher side 
of other measurements of young stars, but consistent to 
within the errors (see following section). The distribu¬ 
tion of metallicities is peaked at ~-0.6dex but with a high 
metallicity tail (see Fig. ED- We also find that the average 
value of £ is slightly lower in the SMC than for th e LMC 
and for the Per OBI stars studied in iGazak et all (2014) 
(£smc = 2.8 ±0.3 km s -1 , as opposed to 3.3 ±0.5 km s -1 
and 3.6 ± 0.5 km s -1 for the LMC and Per OBI respec¬ 
tively, see Fig. E}. Though there does seem to be a 
weak trend of decreasing £ at lower metallicities, this 
is not borne out by a study of RSGs in NGC 6822, 
a galaxy found to have similar metallicity to the SMC 
dPatrick et, al.ll2015lh 

Finally, we note that we do not see any systematic 
trends with the fitted parameters, for example R and 
[Z], which would indicate obvious systematic errors. 

4.1. Comparison of temperatures with SED fits 

In D13 we showed that the measurements of RSG effec¬ 
tive temperatures using the TiO bands were unreliable, 
and presented a more robust way of measuring RSG tem¬ 
peratures using their optical-NIR spectral energy distri¬ 
butions (SEDs). Since the sample of the D13 is identical 
to that used here, it is natural to investigate how the ef¬ 
fective temperatures measured here from the diagnostic 
J-band lines compare to those measured from SED fits. 

In Fig. [5] we plot the T e g measurements of the current 
study against those from the SED fits of D13. The means 
and standard deviations of the quantity AT(= Tj— Tsed) 
are 160 ± 110K and 48 ± 68K for the LMC and SMC 
respectively. This tells us that overall our MARCS mod¬ 
els with non-LTE corrections are able to simultaneously 
match the optical-NIR continuum as well as the strengths 
of the atomic lines in the J-band. Though the agreement 
between the two methods is excellent for the SMC stars, 
there is a systematic offset for the LMC stars, which for 
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Fig. 8. — Comparison of the temperatures derived in this study, T e ff(J), with those from fits to the stellar SEDs (left panel) and from 
fits to the TiO bands in the optical (right panel). 


the two stars with the highest SED temperatures is 3-4cr. 
At the present time we have no explanation for this, or 
whether indeed this offset is significant. A planned study 
of stars at higher and lower metallicities will address this. 

For completeness, in Fig. [5] we also plot the tempera¬ 
tures from the J-band fits against those derived from the 
TiO band strengths, also from D13. We find little cor¬ 
relation between these two, with the TiO temperatures 
being on average 300 K cooler, mirroring the findings in 
D13. 

Since there is a small discrepancy between the J-band 
effective temperatures and those determined from SED 
fits, it is natural to ask what systematic effect this may 
have on the derived metallicities. To answer this ques¬ 
tion we reanalysed the spectra in the same way as in 
Sect. 13.31 but instead of allowing temperature to be a 
free parameter we fixed it to be the same as the SED 
temperatures from D13. The fitted metallicities of each 
star was typically within ler of those quoted in Table [2} 
while the mean metallicities of the LMC and SMC, which 
now were found to be -0.36 and -0.50 respectively, were 
stable to within ±0.03dex. 

5. DISCUSSION: A COMPARISON WITH OTHER 
METALLICITY STUDIES OF THE MAGELLANIC 
CLOUDS 

The average metallicities of the two galaxies are listed 
in Table[2] For the LMC we find [Z] = —0.37 ± 0.14, 
and for the SMC we find [ Z ] = —0.53 ± 0.16. These 
are relative to an a dopted Solar value of Z 0 =O.O12 
(|Grevesse et al.ll2007l) . 

Below we discuss these results in the context of other 
metallicity studies of the Magellanic Clouds. In this 
discussion we restrict ourselves to studies of young 
(^lOOMyr) stars, since these should be the most com¬ 
parable to the results presented here. Compilations of 
these studies are presented in Tables [3] and 0] We also 
note that in our abundance analysis we have kept all met¬ 
als scaled to their Solar values relative to each other. We 
will explore the possibilities of non-Solar a-to-iron ratios 
at the end of this Section. 



Fig. 9. — The best-fit microturbulent velocities as a function of 
metallicity. The average values and t heir standard dev iations for 
the three environments - Per OBI (|Gazak et al.1120141) . and the 
LMC and SMC (this work) - are indicated by the bold symbols 
and thicker error bars. 

5.1. Abundance studies of the LMC 

The most rece nt comparable studie s to that presented 
here a re those of lDufton et a l. (2006) and lTrundle et all 
(2007). These authors studied samples of B stars in clus¬ 
ters and associations in both Magellanic Clouds, using 
non-LTE spectral analysis^- For the LMC, they stud¬ 
ied stars in two clusters, Nil and NGC 2004. In gen¬ 
eral they found abundances of Fe and O which were - 
0.3dex with respect to Solar, and abundances of Si and 
Mg which were lower by ^O.ldex relative to H. A sam¬ 
ple of F su perg iant s were the subject of an LTE study 
by iHill et ah] ( 199511 . finding [0,Mg,Fe] abundances ~- 
0.25dex relative to H, while Si was super-Solar. A re- 

5 The Fe lines were treated in LTE in iTrundle et all (120071 ). 
though non-LTE effects are thought to be small i n the paramete r 
range occupied by the stars in their sample (IThompson et al.120081) . 
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analysis of the same stars by lAndrievskv et al.1 (|200lH 
found an non-LTE-corrected [O] abundance of -0.16dex, 
and [Fe]=-0.34 (LTE). These results are consistent with 
iHill et ahl ( 1995 ). but serve as an illustration of the un¬ 
certainties. Two separate studies of ceph eids fo und abun¬ 
dances in the range - 0.1 to -0.2dex dLuck et al.1 119981 
iRomaniello et al.ll2008l) . 

In general, these previous studies converge at an aver¬ 
age metal abundance [Z] of -0.3±0.1. Our average value 
from our sample of RSGs is therefore perfectly consistent 
with these earlier results. 


5.2. Abundance studies of the SMC 


For this galaxy there is a larger degree of disparity 
between estimates of the metallicity from a variety of 
abundance probes. In Table 0] we summarise selected 
abundance determinations from young stars in the SMC. 
These studies have looked at abundances of individual 
elements rather than the global metallicity Z. We have 
concentrated on the elements O, Si, Mg and Fe (where 
available) since these elements should be relatively un- 
affect ed by st ellar evolution for all but the most evolved 
stars dBrott et al.ll2011l l. 

For studies of hot stars, we consider only those which 
include the effects of non-LTE, which are known to 
be non-negligible in hot star atmospheres. Studies of 
main-sequence B st ars in the yo u ng cl u ster NGC 330 
(iKorn et al.1 l2000t iLennon et al.1 120031: iTrundle et, all 
l2007fl have found O abundances in the range -0.4dex to - 
0.8dex relative to the Solar value of log(Q/H) + 12 = 8.66 
(jGrevesse et alJl2007il . Since one would expect the abun¬ 
dances in a young cluster to be uniform this serves to 
illustrate the uncertainties in this type of work. Abun¬ 
dances of Mg and Si range from -0.7 to -0.9dex com¬ 
pared to Solar. A study of another young cluster, 
NGC 346, and the isolated B star AV 304, found abun¬ 
dances for these three elements that were consistent 
with the NGC 330 studies (-0.5 to -0.8dex for O, Si 


and Mg). Analysis of ev olved hot stars, B supergiants 
(ITrundle fc Lennonll2005l l and A supergiants (|Vennlll995l 
I1999R again revealed similar values, though the latter is 
an LTE study. 

Previous analysis of cool evolved massive stars have 
not included non-LTE effects, but are included here for 
completeness, and are discussed within the context of 
our recent inves tigations into non-LTE corre ctions for 
cool supergiants (IBergemann et al.l[2012L [2013 1. Another 
study of NGC 330, this time of the K supergiants, re¬ 
vealed abundance levels which were slightly lower than 
those obtained from analysis of hot stars, with abun¬ 
dances of O, Si, Mg and Fe all between 0.8-0.9dex below 
Solar (iHilll Il999ll . This could in p a rt be explained by 
non-LTE effects: iBergemann et all (|2012j l showed that 
for certain Fe lines the assumption of LTE would result 
in abundance underestimates > O.ldex at temperatures 
of >4400K, though the effect seeme d to be the opposite 
for Si lines (iBergemann et ahll2013lh The analysis of K 
supergiants in the SMC field were higher, be tween 0.5- 
0.8dex below Solar ([Hill et al.ll 19971 iHillll I997li . 

Finally, from analysis of Cepheids, the abundances 
of O, Si, Mg and Fe were found to be 0.5-0.6dex sub- 
Solar dLuck et al.ll 1998l) . though a more recent study by 
IRomaniello et al.l ( 2008T ) found an iron abundance of - 
0.7dex relative to the Solar value of log(Fe/H) + 12 = 


7.45. Both studies assumed LTE. 

To summarise, young star abundance determinations 
seem to have a somewhat large spread, typically finding 
between 0.5-0.8dex below Solar. In this context, our re¬ 
sult of [ Z] = —0.53 ± 0.16 is at the upper boundary of 
these other results, but consistent within the lcr errors. 

5.3. Departures from Solar-scaled a/Fe 

Throughout this study we have consistently assumed a 
value of [a/Fe]=0.0. Our literature study of analyses of 
young stars (Tables [3] and [4} supports this assumption, 
appearing to show that [a/Fe] for the two Magallanic 
Clouds is Solar to within iO.ldex, with the results of the 
FLAMES survey of massive stars indicating that perhaps 
[Si/Fe] and [Mg/Fe] were depleted by 0.1-0.2dex with 
respect to Solar. However a com pre hens ive study of BA 
supergiants in the SMC bv lSchiller : (j20ld ) concluded that 
[a/Fe] was en hanced by O.ldex with respect to the Solar 
abundances of lGrevesse et al.1 (|2007H . 

Here, we now explore the effect that a small de¬ 
parture from Solar [a/Fe] would have on our results 
if Solar [a/Fe] were incorrectly assumed. To do this, 
we constructed a model spectra with parameters [Z]=- 
0.5, T e ff=3900K, £=3km s -1 , loggi=0.0, and with 
[a/Fe]=±0.2. We then analysed these spectra in the 
same way as in the rest of this work using the models 
with Solar [a/Fe]. 

In Fig. [TUlwe plot the degeneracy between [a/Fe] and 
the model parameters [Z], £ and T e g (there was no de¬ 
tectable degeneracy with log g). In the bottom panel we 
see there is a small trend between T e g and [a/Fe]. This 
can be explained by the fact that the relative strengths 
of species of lines of different excitation potentials are a 
diagnostic of temperature (see Sect. 13.4.1?)) . By altering 
the relative strengths of, for example, Fei and Sil lines, 
we may mimic the spectrum of a star with a slightly dif¬ 
ferent temperature with Solar-scaled abundances. This 
effect however is small increasing [a/Fe] by 0.2dex re¬ 
sults in a change of inferred T e g of only ~50K. 

It is a similar situation for microturbulence. This pa¬ 
rameter is sensitive to the relative strengths of the strong 
and weak lines. In principle this could be estimated from 
the lines of jus t one element, for example Fei as discussed 
in Sect. 13.4.11 In practise, our holistic ^-minimization 
approach considers all lines together to increase the pre¬ 
cision on each parameter, which results in the degener¬ 
acy seen in the middle panel of Fig. [TU] An increase in 
[a/Fe] of 0.2dex would result in f being underestimated 
by ~0.3km s -1 . 

The effect on overall metallicity is predictable ■ alter¬ 
ing the abundances of the a elements, which constitute 
six of the eight diagnostic lines, results in an inferred 
metallicity which is of order that by which the a ele¬ 
ments are altered. The small changes in T e g and £ are 
required to alter the strengths of the Fe I lines relative to 
those of the Si I and Ti I lines in order to provide a better 
fit to the input spectrum. 

The dashed vertical lines in Fig. ITOl illustrate the possi¬ 
ble deviations from Solar [a/Fe] indicated by our litera¬ 
ture search (±0.1dex). If [ a/Fei were to be enhanced by 
O.ldex, as suggested bv lScliilleii (|2010T) . we would expect 
to see overall metallicities enhanced by a similar amount 
with respect to other studies, whilst we might also see 
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[a/ Fe] 

Fig. 10. — The impact on the derived physical properties of a star 
with non-Solar [a/Fe] ratios when its spectrum is analysed with a 
grid of models with Solar [a/Fe]. The panels show the inferred 
metallicity ( top ), microturbulence ( middle ) and temperature (bot¬ 
tom) as a function of the input star’s [a/Fe] ratio, for two global 
metallicities. See text for details. 


average T e g and £ values which were slightly higher if 
the average values of these quantities were independent 
of metallicity, which may or may not be the case. Our 
results for the SMC do show slightly higher abundances 
than other works by around O.ldex, while the average £ 
value does seem to be lower by a few tenths of km s . 
This is circumstantial evidence, albeit rather weak, for a 
small [a/Fe] enhancement in the SMC. 

In the near future we plan to undertake a more thor¬ 
ough analysis of [a/Fe] ratios in these galaxies. This 
will be possible once we have implemented non-LTE cor¬ 
rections to the two Mg I lines in our spectral window 
(Bergemann et al., in press), with the increased number 
of diagnostic lines and a third a-element enabling us to 
separate [a/Fe] from [Fe/H] in our abundance analysis. 

6. SUMMARY & CONCLUSIONS 

We have presented a metallicity study of the Large and 
Small Magellanic Clouds (LMC and SMC respectively) 
using VLT/XSHOOTER near-IR spectroscopy of sam¬ 
ples of Red Supergiants (RSGs). Such stars are young 
and their abundances of Mg, Ti, Si and Fe accurately 
reflect those in the gas phase. We concentrate our anal¬ 
ysis on a narrow window in the J-band, where molecular 
absorption is weak, and where we have shown previously 
that accurate stellar abundances may be obtained. Our 
results can be summarised as follows: 

• Our analysis of the two samples of stars reveal 
metal abundances of elements heavier than He rel¬ 
ative to Solar of [Z]=-0.37±0.14 for the LMC and 
[Z]=-0.53±0.16 for the SMC. Both results are con¬ 
sistent with other studies of young massive stars 
in the literature, though the SMC result is at the 
high end of that found by other comparable works 
(^0.5-0.8dex below Solar). 

• We find best-fitting temperatures which are con¬ 
sistent with those from fits to the optical-infrared 
spectral energy distribution, though there is a small 
systematic offset of marginal statistical significance 
for the LMC stars, 160±110K. 

• The average microturbulent velocities £ in the LMC 
(3.3±0.5km s -1 ) are consistent with those found 
by ourselves in Galactic stars. The average for 
the SMC stars is slightly lower at 2.8±0.3km s -1 . 
Though a low-significance result, we offer two pos¬ 
sible explanations: firstly, it may be indicative of 
the physics of convection at low metallicity; or sec¬ 
ondly, it could be a systematic effect caused by the 
SMC RSGs having slightly super-Solar [a/Fe] ra¬ 
tios. 

In the near future we will explore the effect of non- 
Solar [a/Fe] ratios on our results by incorporating non- 
LTE corrections to the J-band Mg I lines into our analy¬ 
sis. The increased number of diagnostic lines and atomic 
species will then enable the ratio of a-elements to be 
considered separately to the Fe lines. We will also ex¬ 
tend our metallicity baseline to lower [Z] systems such 
as WLM, to further investigate any potential trend of £ 
with decreasing [Z]. 
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TABLE 3 

Average LMC abundances 


Target (s) 

B dwarfs 

Nil -0. 

NGC 2004 -0. 

Cool Supergiants 
F supergiantsf -0. 

” ” (same sample) -0. 
Cepheidsf 
Cepheidsf 


foi 

f®l 

34zb0.08 

-0.38zb0.07 

27zb0.18 

-0.36zb0.03 

22zb0.08 

+0.08zb0.10 

16zb0.08 

- 

-0.01 

-0.12 


fMgl[Ffe] 


-0.47±0.09 -0.22zb0.10f 

-0.45±0.10 -0.29zb0.12f 

-0.32zb0.08 -0.22zb0.06 

-0.34zb0.15f 
-0.14zb0.12 
-0.27zb0.12 


Ref. 


(2006) 

^Tnmdle^^iL (2007) 

j Hin_et r aL (1995) 
Midrievsk^et^aL (200L) 
i Luck = et_aL (1998) 
Romanieno = et_^L (2008) 


Solar 8.66 7.51 7.53 7.45 Grevesse et al. (2007) 

fStudy did not account for non-LTE effects on diagnostic lines. 


TABLE 4 

Average SMC abundances 


Target(s) 

roi 

rsii 

fMg| 

m 

Ref. 

B dwarfs 

NGC 330 

-0.78zb0.20 

-0.67zb0.02 

-0.80zb0.10 


Trundle et al. (2007) 

55 55 

-0.69zb0.13 

-0.91zb0.32 

-0.91zb0.14 

- 

Lennon et al. (2003) 

55 55 

-0.42zb0.30 

-0.68zb0.40 

-0.69zb0.30 

- 

Korn et al. (2000) 

NGC 346 

-0.61&0.10 

-0.71zb0.05 

-0.76zb0.07 

- 

Trundle et al. (20071 

AV 304 

-0.54zb0.13 

-0.76zb0.16 

-0.79zb0.20 

- 

Hunter et al. (2005) 

Hot Supergiants 
field B SGs 

-0.54zb0.13 

-0.75zb0.18 

-0.69zb0.14 


Trundle & Lennon (2005) 

field A SGsf 

-0.52zb0.06 

-0.54zb0.18 

-0.78zb0.10 

-0.80zb0.15 

Venn 11995. 1999) 

Cool Supergiants 

K SGs (NGC 330)f 

-0.9 

-0.8 

-0.9 

-0.8 

Hill 11999) 

field K SGsf 

-0.7 

-0.8 

-0.5 

-0.6 

Hill et al. 11997): Hill (1997) 

Cepheidsf 

-0.6 

-0.5 

-0.6 

-0.6 

Luck et al. (1998) 

Solar 

8.66 

7.51 

7.53 

7.45 

Grevesse et al. (2007) 


fStudy did not account for non-LTE effects on diagnostic lines. 
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